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We present a study of transport in graphene devices on polar insulating substrates using a tight- 
binding model. The mobility is computed using a multiband Boltzmann treatment. We provide 
the scaling of the surface polar phonon contribution to the low-field mobility with carrier density, 
temperature, and distance from the substrate. At high bias, we find that graphene self-heating effect 
is essential to account for the observed saturated current behavior. We predict that by optimizing 
the device cooling, the high bias currents can be significantly enhanced. 



The excellent transport and optical properties of 
graphene [l[ provide strong motivation for research into 
possible applications of this material in nanoscale elec- 
tronics and optoelectronics @, 3 3- The electrostatic 
modulation of the graphene channel through gates yields 
very promising two-dimensional field-effect devices for 
analog and radio- frequency applications Q . Such de- 
vices should ideally be operated in the saturation limit 
Q. Indeed, recently it has been shown that the cur- 
rent saturates as the source-drain field is increased to 
a few Volts per micron @, S, @| ■ While elastic scatter- 
ing determines the rate at which current increases with 
the applied bias, the current saturation process has been 
attributed to either the inelastic scattering of electrons 
by surface polar phonons (SPP) in the polar substrates 
0, [i| or the intrinsic graphene optical phonons [§] . Sig- 
nificant heating of the graphene devices operated under 
high bias conditions has also been observed by Raman 
spectroscopy 0. However, theoretical studies of such 
temperature effects on the inelastic scattering and de- 
vice performance optimization by substrate engineering 
and thermal management are lacking. 

In this Letter, we evaluate diffusive transport prop- 
erties of graphene on SiC>2, H1O2, and SiC polar sub- 
strates by solving the Boltzmann transport equation 
(BTE) in the presence of both intrinsic graphene phonons 
and substrate SPP phonons. The electronic structure of 
graphene is described by a 7r-orbital tight-binding model 
with a hoping parameter to — 3.1 eV, which gives a 
Fermi velocity vf — {V2>/2)toa/h ps 10 6 m/s, where 
a = 0.246 nm is the graphene lattice constant. For 
the electron-phonon scattering we use the Su-Schrieffcr- 
Heeger (SSH) model [HJ for modulation of the 7r-orbital 
overlap t = to — gSRcc with g = 5.3 eV/A as used in 
the mobility calculations in carbon nanotubes [ill ]. The 
parameters for SPP scattering in different substrates are 
given in Table HI 

The low-field mobility in pristine graphene, in the ab- 
sence of charged impurities and defects, is determined 
by scattering from the graphene phonons. The results 
for the low-field mobility are shown in Fig. [T] Within 
the SSH model there are two transverse (TA) and longi- 
tudinal (LA) acoustic phonon branches whose electron- 
phonon couplings can be approximated by the deforma- 
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FIG. 1: (Color online) Temperature dependence of the low- 
field mobility in graphene on different substrates: green - in- 
trinsic, cyan - on SiC, red - on Si02, blue - on Hf02 at carrier 
density (a) n = 10 12 cm" 2 and (b)n=5x 10 12 cm" 2 . The 
solid curves are fits to Eq. ([5]) with three global fit parameters: 
characteristic sound velocity v p h = 17.3 km/s, energy of the 
optical phonon in graphene /kj op = 0.19 eV, and (3 ~ 0.0115 
in Eq. ©. 



tion potential ID, IH : D ac — Sganj (4\/3) , where the re- 
duction factor k — v% a /(y/2(v^ A —v^i )) w as introduced 
in For the valence phonon model [14| used here we 
obtain D ac m 2.6 eV. The TA and LA modes have differ- 
ent angle dependencies for the electron-phonon couplings 
M k , k+q \ 2 = D 2 ac q 2 h/{2M C Lo q )(l ± cos(3(0 fc + k+q )))/2 
12], where 9 k is given by cos#fc = 5k x /\5k\ and sinflfc = 
5k y /\5k\ for the carrier wavevector k in the vicinity of 
the if -point, i.e. k = K + Sk. Thus, the acoustic phonon 
mobility contribution can be approximated by an angle 
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independent coupling and a characteristic sound velocity 
v p h- The low-field mobility in this model has been calcu- 
lated Il5| and in the large temperature limit it was found 
to be UM : 
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where p m is graphene mass density and n is the carrier 
concentration. 

The two optical phonons at T point have couplings 
l M fe,fe'+g| 2 = ^a/(2M c Wr)(l ± ss' cos{0 k - 9 k ^)) for 
LO '(- sign) and TO (+ sign) modes respectively 0,0], 
where D r = Sg/V? « 11.2 cV/A QjJ, s = 1 for elec- 
trons and s = — 1 for holes. The K-point TO phonon 
mode has an electron-phonon coupling twice as large 
[U E3] with the angle dependence given by \M^l +q \' 2 = 
Dlh/{M C 0JK){1 - ss'cos(6 k - 6 k+q )) 0. The effect of 
the optical phonons (both at T and K) on the low-field 
mobility can be calculated according to [l|| : 
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where N op is the Bosc-Einstein occupation number of op- 
tical phonons and D op = 2Z?r = 22.4 eV/A is an effective 
electron-optical phonon coupling [2(| . 



TABLE I: Parameters for the SPP scattering for graphene 
on Si02, Hf02, and SiC substrates. The surface phonon 
(SO) frequencies are obtained from the bulk longitudinal (LO) 



phonons as ujso = u LO ( l+lj^ ) 



i/-> 





Si0 2 [26] 


Hf0 2 [27] 


SiC[28] 


£0 


3.9 


22.0 


9.7 


Ei 


3.36 


6.58 






2.40 


5.03 


6.5 


hujsoi in meV 


58.9 


21.6 




hcoso2 in meV 


156.4 


54.2 


116.0 


F? in meV 


0.237 


0.304 




Fi in meV 


1.612 


0.293 


0.735 



The SPP scattering affects the temperature depen- 
dence of the mobility in graphene 2l|, |22|, |23[ and carbon 
nanotubes |24f on polar substrates. In graphene it is 
given by [22j, 123J : 
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where zq ss 3.5 A is the van der Waals distance between 
the polar substrate and the graphene flake [25| . The mag- 
nitude of the polarization field is given by the Frohlich 

coupling: F 2 = huls °' v f-^rr - r-Vr ) , where fiw so ,„ is 



a surface phonon energy and £o and are the low- and 
high-frequency dielectric constants of the polar substrate. 
The dielectric constant of air is one. When there are sev- 
eral SPP phonon modes with an appreciable coupling, 
then the low- and high- frequency e are understood as 
an intermediate dielectric functions at uji <C loso,v for e o 



and at u>i ^> u>so,l< for £oo [27 1 



We find that the SPP contribution to the low-field mo- 
bility can be approximated as: 
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which is a non-monotonic function of carrier density n. 
Here ko w y/ (2u>spp/vf) 2 + cm, where value of a « 10.5 
was determined from the calculated low-field mobility de- 
pendence on zq (not shown) and parameter (3 ss 0.0115 
is a global fit parameter used in Fig. [TJ N spp ^ is the 
occupation number of the SPP phonons. 

The calculated low-field mobility from the BTE solu- 
tion in Fig. [T] can be well fitted using Matthiessen's rule: 
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where mobility contributions due to the acoustic, optical, 
and SPP phonons are given by Eq. (JTJ) , ((2|), and (|4|) re- 
spectively. Eq. ([5]) describes BTE results remarkably well 
except for the case of Hf02 at high temperatures, where 
both SPP phonons have similar coupling strengths. 

At high bias, the transport is typically described by 
the saturated current model @,[ll|: j = l+^F/j t > where 
(7 = enp is a low-field conductivity and j sat is a satura- 
tion current. In the full saturation regime only carriers 
around Fermi energy Ep in the energy window Ep±Ml/2 
contribute to the current. The saturated current can be 
readily calculated 0,0,11 and for E F > Ml/2: 
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Here is a characteristic frequency of the phonon re- 
sponsible for the current saturation. 

The current densities as a function of electric field 
are shown in Fig. [2] for graphene on different substrates. 
When phonons are kept in thermal equilibrium at T am (, = 
300 K, the current docs not show full saturation for the 
experimentally relevant source-drain fields up to 2 V//im. 
At low densities, the current shows negative differential 
conductance for scattering by intrinsic graphene phonons 
because of the deviation of the bandstructure from the 
linear band dispersion, similar to the effect of the non- 
parabolicity in carbon nanotubes llj. The current at 



high bias (2 V/ pm) in Fig. Q] is proportional to Ep as 
suggested by equation Eq. ([6]). However, the values of 
HQ w 259,238,150 meV for SiC, Si0 2 , and Hf0 2 sub- 
strates correspondingly are significantly larger than the 
SPP phonon energies in Table fl] 
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FIG. 2: (Color online) Current density electric field depen- 
dence in graphene on different polar substrates: green - in- 
trinsic, cyan - on SiC, red - on Si02, blue - on Hf02 at carrier 
density (a) n = 10 12 cm~ 2 , (b) n = 5 x 10 12 cm -2 , and (c) 
n = 10 13 cm~ 2 . The dashed curves show calculations for 
phonons at T am t — 300 K and solid curves show self con- 
sistent calculations with T = T a mb + j(T)F/g for g = 0.278 
kW/(K cm 2 ) 



Recently, it has been shown [8j that graphene un- 
der high bias can experience a significant self-heating 
with graphene phonon temperatures reaching up to 1000 
K. The temperature was found to be proportional to 
the Joule losses T = T am h + jF/g, where the ambient 
temperature was about T am b = 300 K and g = 0.278 
kW/(K cm 2 ) d]. In the presence of the SPP scatter- 
ing, electrons can give their energy directly to the sub- 
strate SPP phonons [29j, which can be heated, in prin- 
ciple, to temperatures higher than that of the graphene 
phonons. In our "self-heating model" we assume that 
SPP and graphene phonons are heated to the same tem- 
perature which is proportional to the Joule losses found 
self-consistcntly, i.e. j = j(T). As a result of self- heating 
the current densities drop by up to a factor of four at 
high biases, especially at high carrier densities, as shown 



in Fig. [2j Moreover, the current shows true saturation at 
experimentally accessible source-drain fields. The high 
bias currents are still proportional to the Ep in the self- 
heating model; however, the current does not extrapo- 
late to zero at low density. The values of Ml as 61, 49, 39 
meV for SiC, SiC-2, and HfC-2 respectively, are extracted 
from the current values at 2V/^tm using Eq. ([6]) and 
are comparable to Ml reported in 0, H|. In the self- 
heating model, the temperature increase contributes to 
the pinch-off effect at low carrier density (see Fig. as 
was recently observed in Ref. At high density, the 
self-heating with intrinsic graphene phonons is predicted 
here to lead to negative differential conductance similar 
to the effect observed in carbon nanotubes 13011. 
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FIG. 3: (Color online) Modeling of the current-voltage char- 
acteristics measured in Ref. and [3] in graphene devices 
on Si02 substrates. The green and red dashed curves, cor- 
respondingly, are BTE solutions for graphene phonons and 
both graphene and SPP phonons scattering at T am b = 300 K 
including Coulomb scattering (see text). The solid green and 
red curves are self-heating model calculations with graphene 
phonon only and both SPP and graphene phonon scattering, 
respectively, (a) Black circles are measurements from Ref. [9| 
(hole branch). The red solid curves used g = 0.455 kW/(K 
cm 2 ) as a fit parameter. The inset shows calculated current 
density at IV / fim field compared to the Isr (black solid cir- 
cles) from Ref. 0. (b) Black squares are from Ref. The 
solid curves are calculated with g = 0.278 kW/(K cm 2 ) @|, 
which are essentially zero parameter fits. 

In the experiment, the low-field mobilities can be sig- 
nificantly effected by the presence of defect and charge 
impurity scattering [3l|, 32 1. As in Ref. Q we include 
charge impurity scattering in the modeling of the exper- 
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imental I-V characteristics in Fig. [3l following [31 L The 
observed mobility of about 1000 cm 2 /Vs in Ref. 8j can 
be reproduced in our calculations by assuming scatter- 
ing with charge impurities of density rii = 4.5 x 10 12 
cm~ 2 and a smaller rii = 3.5 x 10 12 cnr 2 in the presence 
of SPP scattering in SiC>2 substrate. The carrier den- 
sity was fixed by the gate at n ~ 10 13 cm~ 2 in Q and 
in Fig. [3b we show that the calculated current is signifi- 
cantly larger at high biases in the presence of both intrin- 
sic graphene and SPP phonon scattering if the tempera- 
ture is fixed at room temperature. Most importantly, the 
current does not show the saturation that is observed in 
the experiment. On the other hand, using the experimen- 
tally measured temperatures, our self-heating model with 
g = 0.278 kW/(K cm 2 ) and including the SPP scattering 
very nicely reproduces the experiment. At the same time, 
the self-heating model with only the graphene phonons 
active does not show full saturation even at fields up to 2 
V//im and it overestimates the measured current at high 
biases. 

The four probe FV characteristics at low carrier den- 
sity Q were analyzed using defect and charge impurity 
scattering and intrinsic graphene scattering. While the 
overall agreement between the theory and the experi- 
ment was considered satisfactory in [9(, here, we show 
that the agreement can be improved if the data are an- 
alyzed within the self-heating model in the presence of 
SPP scattering from the SiC>2 substrate. To reproduce 
the low-field mobility of about 11000 cm 2 /Vs in Q, we 
use rii = 4.3 x 10 11 cm~ 2 for the Coulomb scattering. A 
similar value was used in The model that includes 
only intrinsic graphene scattering gives a good agreement 
with the measurements [§] at n ~ 1.72 x 10 12 cm -2 in 
Fig. [3^,. However, at lower carrier density the agreement 
is worse as seen from the inset in Fig. [3^, and as it was 
also found by Barrciro et. al. [j| in Fig. 2d. In the pres- 
ence of SPP scattering, we use an impurity concentration 
of rii = 2 x 10 11 cm~ 2 to get the same low-field mobility. 
While the isothermal calculations overestimate the mea- 
sured current at high bias, the self-heating model with 
g = 0.455 kW/(K cm 2 ) reproduces the I-V characteris- 
tics at a fixed density in Fig. [3^ fairly well. Moreover, 
the agreement of the self-heating model with experiment 
holds even at low carrier densities as can be seen in the 
inset of Fig. [3^, which shows the experimentally mea- 
sured current at minimum dl/dV (3 compared to the 
calculated current densities at 1 V//zm. 

In conclusion, our calculations suggest that SPP scat- 
tering is the likely mechanism for the current saturation 
and that the observed full current saturation can only 
be accounted by the self-heating model. Without self- 
heating, the current densities are predicted to be too high 
for either graphene phonon scattering or SPP scattering. 
Therefore, saturated currents can be enhanced if efficient 
device cooling is applied by the appropriate choice of sub- 
strate and the optimization of the graphenc/substratc 



contact thermal resistance. 

We gratefully acknowledge M. Freitag for providing 
data for Fig. [3)3 and I. Meric for helpful discussions. 
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